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Abstract 

The spectra of light-light and heavy-light mesons are calculated within 
the framework of the QCD string model, which is derived from QCD in the 
Wilson loop approach. Special attention is payed to the proper string dy- 
namics that allows us to reproduce the straight-line Regge trajectories with 
the inverse slope being 27rcr for light-light and twice as small for heavy-light 
mesons. We use the model of the rotating QCD string with quarks at the ends 
to calculate the masses of several light-light mesons lying on the lowest Regge 
trajectories and compare them with the experimental data as well as with 
the predictions of other models. The masses of several low-lying orbitally and 
radially excited heavy-light states in the D, Dg, B, and Bg meson spectra 
are calculated in the einbein (auxiliary) field approach, which has proven to 
be rather accurate in various calculations for relativistic systems. The results 
for the spectra are compared with the experimental and recent lattice data. 
It is demonstrated that an account of the proper string dynamics encoded in 
the so-called string correction to the interquark interaction leads to an extra 
negative contribution to the masses of orbitally excited states that resolves 
the problem of the identification of the D(2637) state recently claimed by the 
DELPHI Collaboration. For the heavy-light system we extract the constants 



*yulia@heron. itep.ru ** nefediev@cfif.ist.utl.pt *** simonov@heron.itep.ru 



1 



A, Ai, and A2 used in Heavy Quark Effective Theory (HQET) and find good 
agreement with the results of other approaches. 

PACS: 12.38.Aw, 12.39.Hg, 12.39.Ki 

I. INTRODUCTION 

The description of the mass spectrum of hadrons is one of the fundamental problems of 
strong interactions. It has been attacked in a sequence of approaches motivated by QCD, 
but still attracts considerable attention. One of the most intriguing phenomena — namely, 
the formation of an extended object, the QCD string, between the colour constituents inside 
hadrons, — plays a crucial role in understanding their properties. In the present paper 
this role is exemplified by spectra of the mass of light-light and heavy-light mesons. In the 
former case we study the role played by the QCD string in formation of the straight-line 
Regge trajectories and discuss the form of the interquark interaction inside light hadrons. 
For heavy-light mesons we find the masses of several low- lying states in the D, Dg, B, and 
Bg meson spectra including orbitally and radially excited ones. 

We calculate and discuss the spin-spin and spin-orbit splittings and compare them to 
experimental and recent lattice data. Special attention is payed to the role of the proper 
string dynamics in establishing the correct slope of the Regge trajectories for both light-light 
and heavy-light states, as opposed to those following from relativistic equations with local 
potentials. 

We remind the reader then that an extra piece of the effective interquark potential, the 
string correction, which is entirely due to the string-type interaction in QCD gives a 
negative contribution to the masses of orbitally excited states. The latter observation allows 
us to resolve the "mystery" of an extremely narrow 1^(2637) state (and a similar one in the 
5-mesonic spectrum) ||^ recently claimed by the DELPHI Collaboration We present a 
reasonable fit for the several lowest states in D- and 5-mesonic spectra using the standard 
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values for the string tension, the strong couphng constant, and the current quark masses. 
We also find the correspondence between our model and Heavy Quark Effective Theory, 
extracting the constants used in the latter approach in the expansion of a heavy-light meson 
mass in the inverse powers of the heavy quark mass. We find analytical formulae for these 
constants and compare their numerical estimates with the predictions of other models. 

The two main approaches used in the numerical calculations are the quasiclassical method 
of solving the eigenenergies problem and the variational one based on the einbein field 
formalism. The accuracy of both methods is tested using exactly solvable equations and 
found to be about 7% at worst even for the lowest states. Possible improvements of the 
method are outlined and discussed. 

The paper is organized as follows. In Section || we give a brief insight into various aspects 



of the einbein field formalism. In Section III the exact spectra of relativistic equations are 



compared to the results of approximate calculations using the quasiclassical and variational 
einbein field methods, as well as the combined one. In Section ^ we discuss the problem 
of the Regge trajectory slopes as they appear from the relativistic equations with local 
potentials and from the string-like picture of confinement. Derivation of the Hamiltonian 
for the spinless quark-antiquark system as well as of the spin-dependent corrections to it 
is the subject of Section 0. The spectra of light-light and heavy-light mesonic states are 
calculated and discussed in Sections and [VI I| , respectively. Section [VII1| contains our 



conclusions and outlook. 

II. EINBEIN FIELD FORMALISM 

In this section we give a short introduction to the method of einbein fields and its 
possible applications to relativistic systems. The interested reader can find a more detailed 
information in HI and references therein. 
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A. Reparametrization invariance and constrained systems 



Historically the einbein field formalism was introduced in [0 to treat the kinematics of 
the relativistic spinless particles. Later it was generalized to the case of spinning particles 
1^ and strings ||^. So the action of a free relativistic particle can be rewritten asQ 

S = /J L(r), L = -mv^ ^ " ^' (1) 

where the dot denotes derivative with respect to the proper time r, /i being the einbein 
field0. The original form of the action can be easily restored after solving the Euler-Lagrange 
equation of motion for the einbein fi which amounts to taking the extremum in the latter. 
Note that the invariance of the initial action with respect to the change of the proper time, 

r^fir), ^>0, /(r,) = r„ f{Tf)=rj, (2) 

is preserved if an appropriate rescaling is prescribed to fi: 

- f^/f- (3) 

The latter invariance means that one deals with a constrained system. For the free 
particle the only constraint defines the mass shell, 

-m^ = 0, (4) 

or, in presence of the einbein field /z, 

.^0. H-El^, (5, 



^In the path integral formalism this transformation is based on the following relation 
J Dfi{T) exp ^- J dr ^ ^ + ^) ) ~ ^xp ^- J drVab] . 

^Usually e = ^ is referred to as the einbein Q. 
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with TT being the momentum canonically conjugated to /i and H being the Hamiltonian func- 
tion of the system (in case (^) it identically vanishes). The requirement that the constraint 
TT = be preserved in time returns one to the mass-shell condition (^): 

= vr = {vri/} = — = ^-^~p -m. (6) 

To make things simpler, one can fix the gauge-like freedom (j^) identifying the proper 
time T with one of the physical coordinates of the particle. The most popular choices are 

• the laboratory gauge (r = Xq); 

• the proper time gauge (r = (nx) , = with being the total momentum of the 
system) 0; 



• the light-cone gauge (r = |(xo + 2:3) = x+), 

which lead to quantization of the system on different hypersurfaces. 

With the laboratory gauge fixed the Lagrangian function ([^) becomes 

.2 ^2 



L = -^-^ + ^, (7) 

so that the corresponding Hamiltonian function reads 

i/ = ^±^ + ^, (8) 

2/x 2' ^ ^ 

and after taking the extremum in /i one ends with the standard relativistic expression 

H = ^f + m'^. (9) 



B. Einbeins as variational parameters 

In the simple example considered above neither the Lagrange nor the Hamilton functions 
of the system contained /i, which allowed one to get rid of /i at any stage by taking the 
extremum in the latter. It is not so for more complicated systems when a change of variables 
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is to be performed which touches upon the einbeins. The velocity corresponding to the 
original degrees of freedom of the system may mix in a very tangled way with those for 



einbeins, so that it is not a simple task anymore to follow the lines d la Dirac |11| to resolve 



the set of constraints and to get rid of nonphysical degrees of freedom. See, e.g., ||g,§,|T^] for 
several examples when such a resolution can be done explicitly. 

Luckily another approach to einbeins is known p|,p!3|. They can be treated as variational 
parameters. Thus one replaces the dynamical function of time /u(t) by the parameter /ig 
independent of r. The eigenstate problem is solved then, keeping /io constant, so that one 
has the spectrum M{„}(/io), where {n} denotes the full set of quantum numbers. Then one 
is to minimize each eigenenergy independently with respect to fi(Q. 

aM|„}(/io) 



= 0, M|„| = M|„| (10) 



Such an approach has a number of advantages. First, it allows one to avoid the tedious 
algebra of commuting constraints with one another following the standard Dirac technique 



|rT|. Second, it allows one a very simple and physically transparent interpretation of ein- 
beins. Indeed, in formulae (|lD and (|^) the einbein fi can be treated as an effective mass of 
the particle; the dynamics of the system remains essentially relativistic, though being non- 
relativistic in form. If m is the current quark mass, then /i can be viewed as its constituent 
mass celebrated in hadronic phenomenology. What is more, the current mass can be even 
put to zero, whereas the Lagrangian approach remains valid in the presence of the einbeins 
and the standard Hamiltonian technique can be developed then. The latter observation is 
intensively used in analytic QCD calculations for glue describing gluonic degrees of freedom 
in glueballs and hybrids |T^|T^. 



^Note that solutions for of both signs appear, but only one of them (hq > 0) is finally left. 
Neglecting the negative solution is the general lack of the einbein field approach and this leads to 
the fact that quark Zitterbewegung is not taken into account (see also the discussion in Subsection 

niDl) . 
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An obvious disadvantage of the variational approach to the einbein fields is some loss of 
accuracy. As a variational method it provides only an approximate solution, giving no hint 
as to how to estimate the ultimate accuracy of the results. Thus in the next section we test 
this method, comparing its predictions with exact solutions of some relativistic equations. 
We consider the accuracy, found to be about 7% at worst, quite reasonable, which justifies 
our consequent attack on the light-light and heavy-light mesons spectra using this formalism. 



III. TESTING THE METHOD 

A. Quasiclassics for the spinless Salpeter equation 

We start from the Salpeter equation for the quark-antiquark system with equal masses 
and restrict ourselves to the zero-angular-momentum case for simplicity: 

{2^pl + m^ + ar^ = M^V'n, (11) 

where the subscript (//) stands for the light-light system. 
The quasiclassical quantization condition looks like 

p^{r)dr = Ti[n + -y n = 0,l,2,..., r+ = , (12) 

where the integral on the l.h.s. can be worked out analytically, yielding 



^(Mr)'-4m2 + Mf) . 3. 



Mf )^(m(")) - 4m2 - Am" In ^ — = Aai: [n + -V (13) 

or approximately (m <^ ^/a) one has 

(m(")) ^ = 4vra + I) + 2m^ In ![^(!i±M + . . . . (14) 



Solution (|l^) becomes exact in the limit m = 0, whereas for a nonzero mass the leading 
correction to the linear regime (M^^^)^ ~ n behaves like 

^.of^m^^U!^. (16) 



For a heavy-light system one has the Salpeter equation 

{JVI + ^^ + = ^)Vn, (16) 

where M'^^'^ denotes the excess over the heavy particle mass M. Similarly to (0) one finds 
then 



Mf ')^(m('^^)) -m^-m' In ^ ^ — 



2aix{n+-\, (17) 



m V 4 

and formula (pIS]) holds true in this well. 

Comparing the results of the WKB method with the exact solutions of equation ([111 
(rows M„(WKB) and M„ (exact) in Table |), one can see that the error does not exceed 3 



4% even for the ground state. See also ||T6[ where the WKB method is tested for light-light 



mesons. 



B. Quasiclassics for the one-particle Dirac equation 

As a next example we discuss the one-particle Dirac equation with linearly rising confin- 
ing potential |1T7| : 

{ap + Pirn + U) + V)^n = Sn^n. (18) 



The WKB method applied to this equation gives p8| . [T9[] 

Kw\ ( 1 

— dr = TT n + - 



— )rfr = 7rfn + i), n = 0,1,2,..., (19) 



where 



P = y (e - V^)2 - ^ - (m + f/)2, (20) 
1 1 U' -Y 

w = — 



2r 2m + U + e -V 

, , 1 
11^2 
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For the most interesting case of purely scalar confinement {V = 0, U = err) an approxi- 
mate quasiclassical solution was found in |T9[ (m = 0): 



A detailed comparison of the results of the WKB method and those following from the 
recursive formula ( pT]) with exact numerical solutions to equation (|I8|) is given in |^ . Here 



we only note that the coincidence of the three numbers is impressive as even for the lowest 
states the discrepancy does not exceed 1%. 

C. Quasiclassical variational einbein field (combined) method for the spinless 

Salpeter equation 

Finally we combine the two methods discussed above and apply the WKB approximation 
to the Hamiltonian of a relativistic system with einbeins introduced as variational parame- 
ters. Then the resulting quasiclassical spectrum is minimized with respect to the einbeins. 
Thus we have a powerful method of solving the eigenvalues problem for various relativistic 
systems which we call "combined." Let us test the accuracy of this method first. 

We start from the Salpeter equation ([TT| ) for the light-light system and introduce the 
parameter /io as described in Section |T|: 

I T)^ -)- 777,2 

Hi = 2Jpj + m^ + ar — > = + /^o + c^r. (22) 



In what follows we consider the massless case substituting m = into (^21). 

We give the analytic formulae for the spectrum of the Salpeter equation ([TT|) obtained 
using the quasiclassical approximation for the Hamiltonian Hi (following from equation 



(p!3|) for m = 0), the exact solution for the Hamiltonian H2 minimized with respect to the 
einbein field and the result of the combined method when the Bohr-Sommerfeld quantization 
condition is applied to the Hamiltonian H2 and the ultimate spectrum is also minimized with 
respect to [Iq. 
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M^(WKB) = 47ra(n + ^), (23) 



M^(einbein) = 16a ( -^I^] ' , (24) 



(combined) = + ^) , (25) 

where Cn+i is the (?t.+ l)th zero of the Airy function Aii^z) and counting of zeros starts from 
unity. The extremal values of the einbein field in the latter two cases read 

3/4 



/Zq (einbein) = \fo [ — ) , /Iq (combined) 



3 ; • 2V3 • 

i.e., the effective quark mass is /ig ~ and it appears entirely due to the interquark 
interaction. 

In Table | we compare the results of the above three approximate methods of solving the 
eigenvalues problem for equation (|IT]) with the exact solution. In the last row we give the 
accuracy of the combined method vs the exact solution. Two conclusions can be deduced 
from Table |I[ The first one is that the accuracy of all approximate methods is high enough, 
including the combined method, which is of most interest for us in view of its consequent 
applications to the QCD string with quarks at the ends. The other conclusion is that the 
variational einbein field method gives a systematic overestimation for the excited states 
which is of order 5-7%. 



D. Discussion 



Here we would like to make a couple of concluding comments concerning the numerical 
methods tested in this section, their accuracy and possible ways of their improvement. As 
stated above the combined quasiclassical variational method is of most interest for us, so we 
shall concentrate basically on it. The following two remarks are in order here. 

From Table |I| one can see that the relative error is practically constant, tending to 
the value of 7% for large n. The reason for such a behaviour will become clear if one 
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compares formulae (p3D and (^). Both relations reproduce the same dependence on the 
radial quantum number n, whereas the difference comes from different slopes, Ana in ( p3D 
vs in (pSl). Then for highly excited states the error is practically independent of n and 
can be estimated as 



^ _ M„(combined) - M„(exact) _ M„(combined) - M„(WKB) _ ^ / ^ ~ o 07 (27) 
~ M„(combined) M„ (combined) ~ ~ V ^ ' ^ 

i.e., the ultimate accuracy of the quasiclassical variational einbein field method (combined 
method) appears to be about 7%. Introducing, say, a correcting factor in (^) one could 
overcome the systematic overestimation and reproduce the spectrum with a better accuracy. 
We shall return to this observation later on when discussing the spectrum of the heavy-light 
mesons. 

Another source of error in the einbein field approach is neglecting the quark Zitterbe- 
wegung (see the footnote on page H). As stated above we neglect the negative sign solution 
for the einbein field fiQ expecting its small influence on the spectrum. Let us give some 
reasoning to justify this action. 

It was demonstrated numerically in ^^ that the contribution of the quark Z-graphs 



into is nearly constant for large excitation numbers and is of order 10%, so that the 
corresponding shift of M behaves like 

AM ~ ^ ~ ^. (28) 



so it is somewhat suppressed. 

Besides, the good agreement of our numerical results with those provided by the lattice 
data and taken from the Particle Data Group can also serve as an a posteriori justification 
of such a neglect. Still some improvements for the einbein field approach are needed to take 
this effect into account. 
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IV. REGGE TRAJECTORIES FOR RELATIVISTIC EQUATIONS WITH LOCAL 

POTENTIALS 



It was observed long ago that the mesonic Regge trajectories are almost linear if the 
total momentum or the radial quantum number is plotted vs the mesonic mass squared ||22|| : 

M^{n,J) = Cnn + cjJ + AM'^, (29) 

where c„ and cj are the (inverse) slopes while AM^ denotes corrections to the leading linear 
regime which come from the self-energy, spin splittings, etc. 

Relations like (^) naturally appear in most of models for confinement, though the (in- 
verse) slopes Cn and cj are different for different models. 

For example the Salpeter equation for the heavy-light system, 

{\/Pr+Pl^ + m^ + i)nl = Mnli)nh (30) 

gives 

cf^\Salpeter) = 4a, cl^^\Salpeter) = 2a, (31) 

where the total momentum J coincides with the orbital one /. 

For the light-light system one easily finds, from (0) by a trivial parameters rescaling, 

cf {Salpeter) = 8a, c^^^\Salpeter) = 4a. (32) 

The one-particle Dirac equation (|18D yields different slopes for different natures of the 
confining force. Thus for the purely vector confinement (potential added to the energy 
termQ) one finds 

c^^UDtrac) = Aa, (33) 
whereas, for purely scalar confinement (potential added to the mass term), 

^We disregard the problem of the Klein paradox here. 
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c?2,i(Dzrac) = 2a. (34) 



In the meantime the spectrum ( pQ]) is expected to follow from a string-like picture of 
confinement which predicts the (inverse) Regge slopes to be 

c'j\string) = 27rcr, c^j^\string) = na. (35) 

One can easily see that none of the relativistic equations considered before gives the 



correct result (^); moreover the discrepancies are rather large (of order 25%). See also ||23 
for a discussion of various models of confinement and the corresponding Regge trajectory 
slopes. 

The reason why relativistic Salpeter and Dirac equations fail to reproduce the correct 
string slope of the Regge trajectories is obvious and quite physically transparent. Indeed, all 
relativistic equations with local potentials have only a trivial dependence of the interquark 
interaction on the angular momentum which comes entirely from the quark kinetic energy. 
Meanwhile, QCD is believed to lead to a string-type interaction between the colour con- 
stituents inside hadrons, whereas the QCD string developed between quarks possesses its 
own inertia and thus it should also contribute to the J-dependent part of the interaction. 
It is this extra purely string-type piece of the interquark interaction to give an extra con- 
tribution to the Regge trajectory slope and to bring it into the correct form of (|35|). This 
statement is proved explicitly in the next section, whereas the string dynamics footprint in 
the heavy-light mesons spectrum is discussed in detail in Section |VI1| . 

V. HAMILTONIAN OF THE QQ MESON 

A. Quark- ant iquark Green's function 

We start from the Euclidean Green's function of the qq pair in the confining vacuum 

G,^ = {¥J){x,y\Ay'^^k^,y\A)),^A, (36) 
where the initial and the final mesonic wave functions 
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^'i^fix^ylA) = ^g{xMx,y)T'~''f^^^{y) (37) 

are gauge invariant due to the standard path-ordered parallel transporter 

y) = P exp (^ig dz^A^^ , (38) 

denote the matrices which might be inserted into the initial and final meson-quark- 
antiquark vertices. 

Integrating out the quark fields in ([36|), one finds, for the mesonic Green's function, 



Ggg = {TrT^f^S,ix,x\AMx,y)T^'^S,iy,y\A)^y,x))A, (39) 

where the trace stands for both colour and spinor indices. We have neglected here the 
1/A^c-suppressed quark determinant, describing sea quark pairs. 

To proceed further we employ the Feynman-Schwinger representation for the one-fermion 
propagators in the external field, fix the laboratory gauge for both particles, 

Xw=ti, X20 = t2, (40) 

and introduce the einbein fields /ii and fi2 by means of the following change of variables (see 



2|,|| for details): 



T 

l^iiU) = T^iioiU), dsiDxio -> Djiiiti), i = 1,2, (41) 
where si,2 are the Schwinger times, T = ^{xq + xq — yQ — yo). 



Then the resulting expression for the mesonic Green's function reads |24 
Ggg = J Dfii{ti)Dfi2{t2)DxiDx2e-'''-'''Tr 



r(-^) (mi - Z))rW (ma - D)x (42) 



\J0 2Hi{ti) ^ ldS^u{Xl{tl)) J \ Jo 2/i2(t2) idSi,u{X2{t2)) J 

with Ki being the kinetic energies of the quarks 
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<^iiu = jiil^ilu — luln)i and S/Ss^i, denotes the derivative with respect to the element of the 
area S. We have also used the minimal area law asymptotic for an isolated Wilson loop, 

(rrPexp {ig dz^A^ ~ exp (-(xS'min), (44) 

which is usually assumed for the stochastic QCD vacuum (see, e.g., ||2^) and found on the 
lattice. Here ^min is the area of the minimal surface swept by the quark and antiquark 
trajectories. 

Looking at (^2]) one can easily recognize the following three main ingredients: the con- 
tribution of the quark, the one of the antiquark, and finally the confining interaction given 
by the string with tension a. One can write, for the latter, 

^min = dt £ dl3^{ww'Y - wV, (45) 

with w^{t,(3) being the string profile function chosen in linear form, 

w^(t, /3) = /Sxi^it) + (1 - (3)x2^, (46) 

thus describing the straight-line string which is a reasonable approximation for the minimal 
surface [0]. 

Finally, synchronizing the quark and the antiquark times (ti = t2 = t) one finds from 
(^) that in the spinless approximation the quark- antiquark meson can be described by the 
Lagrangian 

(47) 

where f = Xi — X2 and n = f/r. Expansion of the surface-ordered exponents in (^) gives a 
set of spin-dependent corrections to the leading regime (|47|) . 

B. Hamiltonian for spinless quarks 



Starting from the Lagrangian (H^ and introducing an extra einbein field z/(t, (3) contin- 
uously depending on the internal string coordinate j3 one can get rid of the square root in 
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(^) arriving at the Hamiltonian of the qq system in the centre-of-mass frame in the form 



1=1 



2j Jo •\2u 2 J 2r2[/ii(l-Cf + /i2C' + /JW/5-C) 



2i ' 



^ ^ /ii + Jo d(3v(3 ^^g^ 

+ ^^2 + Jo d(3v 

Similarly to /i's which have the meaning of the constituent quark masses, the einbein v 
can be viewed as the density of the string energy. In the simplest case of / = one easily 
finds, for the extremal value of 

i^o = (^r; (49) 

i.e., the energy distribution is uniform and the resulting inter quark interaction is just the 
linearly rising potential ar. In the meantime, if / 7^ 0, then the two contributions can be 
identified in the last /-dependent term in (^Sf). Roughly speaking the first two yU-dependent 
terms in the denominator come from the quark kinetic energy. The last term containing 
the integral over /3 is nothing but the extra inertia of the string discussed before. Rotating 
string also contributes to the interquark interaction making it essentially nonlocal, so that 
the very notion of the interquark potential is not applicable to the system anymore. 

Note that the Hamiltonian (|48|) has the form of sum of the "kinetic" and the "potential" 
parts, but this is somewhat misleading, as extrema in all three einbeins are understood, 
so that the ultimate form of the Hamiltonian would be extremely complicated and hardly 
available for further analytical studies. 

Expression ( ^H]) can be simplified if one expands the Hamiltonian in powers of ^/o^/ fi. 
One finds then HI 

H = Ho + Vstring, (50) 

H„ = Hl±^^^A^.r, (51) 



i=l 



2fii 2 
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where Vstring is known as the string correction and this is the term totally missing in the 
relativistic equations with local potentials. Indeed, the Salpeter equation with the linearly 
rising potential is readily reproduced from (|5ll) if extrema in /xi 2 are taken explicitly, whereas 
the string correction is lost. Meanwhile, its sign is negative so that the contribution of the 
string lowers the energy of the system, thus giving a negative contribution to the masses of 
orbitally excited states, leaving those with / = intact. In Section ^ we shall demonstrate 
how a proper account of the string dynamics in the full Hamiltonian ( ^HD brings the Regge 
trajectory slope to the correct value (pSf), whereas in Section |VI1] the string correction ( |52D 
will be demonstrated to solve the problem of the identification of the resonance D[2Q37) 
recently claimed by the DELPHI Collaboration [^0. 



C. Spin-dependent corrections 



Let us return to the quark- ant iquark Green's function (|4^) and extract the nonpertur- 
bative spin-orbit interaction. Following [0,^] one finds 



It follows from |2^j2^ that all potentials Vi{r) (in the notation of contain both 
perturbative and nonperturbative pieces given there in explicit form. One can argue that at 
large distances only the piece (^3]) is left whereas for hght quarks all nonperturbative ones 
may be important (see PP|). 

Now, to have a full picture of the interquark interaction one is to supply the purely non- 
perturbative string-type interaction described by the Hamiltonian (^Sp by the perturbative 
gluon exchange adding the colour Coulomb potential to the Hamiltonian Hq from (|5ll) and 
calculating the corresponding spin-dependent perturbative terms in addition to the potential 
(0). The result reads 
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where we have also added the overall constant shift Cq and 

T/ ^^'^ c M /^nM2 (SiL S2L\ k ( 1 1 \ SiL K ( 1 1\ S2L 
3/iiyU2 2r \^ /if /i^ y H \2/ii /i2 / yUi H \2^2 /ii / 

2 

+ (3(5irr)(S2n) - (5i52)) + (SL) (2 - ln(/ir) - 7^;), 7^ = 0.57, (55) 

with K = Icto, u = We have also added the term of order which comes from 

one-loop calculations and is intensively discussed in the literature ||30| , |31|j29[| . It is important 
to stress that Cq is due to the nonperturbative self-energy of light quarks, which explains 
the later numerical inputs. 

An important comment concerning the expansion ( |55D is in order. Up to the last term 



the expression (^) coincides in form with the Eichten-Feiberg-Gromes result [Q, but we 
have effective quark masses /ij in the denominators instead of the current ones rrii. Once 



/^j ~ Y ii^) + rnf > rrii or even fii rrii, then the result (|55|) is applicable to the case of 
light quark flavours, when the expansion of the interaction in the inverse powers of the quark 
mass rrii obviously fails. The values of yu's are defined dynamically and differ from state to 
state (see |2^ for details). 

The Hamiltonian (|5^) with spin-dependent terms ( ^5]) will be used for explicit calcula- 
tions for heavy-light mesons. In the case of light-light states one should include additional 



nonperturbative spin-dependent terms (see p6[ and references therein). The masses of light- 
light mesons listed in Table |I| have been calculated from the Regge trajectories which do 
not take into account spin-dependent terms and we give them for the sake of comparison. A 
more detailed calculation for the light mesons taking these effects into account can be found 
2y|. 



m 



VI. SPECTRUM OF LIGHT-LIGHT MESONS 
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A. Angular-momentum-dependent potential and Regge trajectories 



Starting from the Hamiltonian ( ^HD we stick with the case of equal quark masses m 



h=pi±^+,+ — — +^ r'^dp. (56) 



The extremal value of the einbein field v can be found explicitly and reads [20 



u,{(3) = ^ , (57) 



2 

where y is the solution of the transcendental equation 

= -^(arcsin y - yJl-y^) + — , (58) 
ar ^ ar 



and L2 = /(/ + !). 

For large angular momenta the contribution of the quarks (the last term on the r.h.s. 
of (|58D ) is negligible so that the maximal possible value y = 1 is reached, thus yielding the 
solution for the free open string [Q] (see also the second entry in |]T2|). 

With the extremal value z/q from (|57|) inserted, the Hamiltonian (|56D takes the form 



pi + rm?' . . ar . , , o ,rr^\ 

n = --^ h /i(r) H arcsm y + MTjy , (59) 

y 



with y defined by equation (0). The last two terms on the r.h.s. of equation (^Uj) can be 
considered as an effective "potential" 



(TV 

U{fi,r) = — arcsin y + fx{T)y'^, (60) 



which is nontrivially / dependent. In Fig.l we give the form of the effective potential ( |60D 
for a couple of low-lying states (solid line). It has the same asymptotic as the naive sum 
of the linearly rising potential and the centrifugal barrier coming from the kinetic energy of 
the quarks (dotted line). In the meantime it differs from the latter at finite distances. The 
only exception is the case of zero angular momentum which should be treated separately 
and leads to the linearly rising potential for any interquark separation. 
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B. Numerical results 



Following the variational einbein field method described and tested above, we start from 



the Hamiltonian (59) and change the einbein field /j, for the variational parameter /iq [20 
so that one has 

H = ^^^^ + fio + U{fio,r), (61) 

U{fio,r) = — arcsin y + fioy'^. (62) 
Then the quasiclassical method applied to the Hamiltonian (^) gives 



r+ 



Pr{r)dr = TT (^n + , (63) 



with 

Pr-(r) = v//io(M-/io-f/(/xo,r))-m2. (64) 

The eigenvalues M„;(/io) for m = were found numerically from (|63|), (Q) and the 
minimization procedure with respect to /iq was used then. Results for M„; are given in 
Table ^ and depicted in Fig. 2 demonstrating very nearly straight lines with approximately 
string slope (27rcr)^^ in / and twice as small slope in n. Note that it is the region of 
intermediate values of r that plays the crucial role in the Bohr-Sommerfeld integral (|63|), 
i.e., the region where the nontrivial dependence of the effective potential U{fiQ,r) on the 
angular momentum is most important (see Fig.l). 

In Table |ITl| we give a comparison of the masses of several light-light mesonic states 
extracted by means of the numerical results from Table |I| with the experimental data and 



theoretical predictions taken from We have fitted our results to the experimental 



spectrum using the negative constant AM^ (see equation (p9D). 
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C. Discussion 



Let us recall the results obtained for the light-light mesons and discuss problems con- 
nected to the given approach. The net result of the current section is the /-dependent 
effective interquark potential which gives the naive linearly rising interaction only for I = 0. 
It was observed long ago P5|,3 that for large angular momenta the quark dynamics is neg- 
ligible and the slope (P3D naturally appears from the picture of open rotating string. In the 
present paper we find that for massless quarks even the low-lying mesonic states demonstrate 
nearly straight-line Regge trajectories with string slope (PB]). 

One problem clearly seen from Figs. 2, 3 is the leading trajectory intercept /q = I (M^ = 0). 
To reproduce the experimental intercept around -0.5 (see Fig. 3) starting from the theoretical 
one -1-0.5 (see left plot in Fig.2) one needs a large negative constant added either to the 
Hamiltonian (|48D (see, e.g., Cq in equation (^If)) or in the form of AM^ directly in (|9D 
(see also Table |ITTD . Once the first way might violate the linearity of the Regge trajectories, 
then one should expect QCD to prefer the second one, though the first way remains more 
attractive from the practical point of view and will be used in calculations of the heavy-light 
mesons spectrum in the next section. 

Another problem is that one of the most intriguing questions of mesonic spectroscopy, the 
7T—p splitting (and a similar problem in the strange sector) cannot be addressed in our model. 



Taking the exact solution of the spinless Salpeter equation (|TT]) with n = / = (see Table | 
with an appropriate rescaling from cr = 0.2GeV^ to a = O.nGeV"^) one finds for the p mass 
squared the value of order l.TGeV"^ which does not violate the linearity of the trajectory 



(see the circled dot in Fig.2). If the overall negative shift with wjAM^I = 1126MeV"^ (see 



Table p4 ) is applied to this state, then one arrives at a p-meson mass about 775MeV, i.e., a 
value very close to the experimental one. Note that we have practically coinciding constants 
for the p- and a-meson trajectories (see the caption for Table [ITD , which supports the idea 
that AM^ can be associated with quark self-energies. 

Meanwhile, one cannot pretend to describe pions (kaons) in the same framework as their 
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Goldstone nature is not implemented in the current model. In realistic quantum-field-theory- 
based models each mesonic state possesses two wave functions which describe the motion 



forward and backward in time of the qq pair inside the meson The backward motion is 
suppressed if at least one of the quarks is heavy, for highly excited states and in the infinite- 
momentum frame. For the chiral pion, which is expected to be strictly massless in the chiral 
limit, the two wave functions are of the same order of magnitude (see, e.g., [^] for an explicit 
pionic solution in QCD2), so that none of them can be neglected. This explains why the 
naive estimate for the pion mass lies much higher than the experimentally observed value of 
lAOMeV. For the first excited state, p meson, this effect is already suppressed, though one 
still has to be careful to neglect the backward motion of the quarks. The progress in this 



direction was achieved in recent papers by one of the authors (Yu.S.) [T^, where a Dirac- 
type equation was derived for the heavy-light system and the properties of its solutions 
were investigated. This new formalism is expected to allow consideration of pionic Regge 
trajectories as it has the chiral symmetry breaking built in. 

VII. SPECTRUM OF HEAVY-LIGHT MESONS 



All results obtained for the light-light mesons in Section ^ can be reproduced for the 
heavy-light states, so that in the one-body limit the Regge trajectories with the correct 
string (inverse) slope na are readily reproduced. Meanwhile, the aim of this study is to take 
into account corrections to the leading regime which come from the spin-dependent terms 
in the interquark interaction as well as those due to the finitness of the heavy quark mass. 
Corrections of both types are important for establishing the correct spectra of D and B 
mesons which are the main target of the present investigation. 

A. Spectrum of the spinless heavy-light system 



In this subsection we study the spectrum of the heavy-light mesons, disregarding the 
quark spins. This amounts to solving the Schrodinger-like equation for the Hamiltonian Hq 
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from (|5J). Note that to this end one needs to know the nonrelativistic spectrum in the 



potential which is the sum of the hnearly rising and Coulomb parts [26 



ci2 



+ \x\ 



X 



\x\ 



(65) 



where 



2/3 



4 /ilyU2 
6 /Ui + /i2 



If solutions of (1^) for x\ are known as functions of the reduced Coulomb 

potential strength A, then one can find the following expressions for the extremal values of 
the einbeins (constituent quark masses): 



f,,{X) = ^ml + A'{X), /i2(A) = V^i + A2(A), f^{X) = ^-V^ (^] 



3/2 



(66) 



with A(A) given by 



A^(A) 



aX 

3k, 



a + 2A 



da 



dX 



The definition of the reduced einbein field fi via fii and fi2 leads to the equation defining A 

/ii(A)/i2(A) 



/x(A) 



(67) 



/ii(A)+/i2(A)' 

Technically this means that one should generate self-consistent solutions to equations 
(|65|) and (^) which are subject to numerical calculations In Table |^ we give such 

solutions for several radial and orbital excitations in D—, Dg—, B—, and 5^— mesonic spec- 
tra. We use the standard values for the string tension, the strong coupling constant, and 



the current quarks masses. Note that is chosen close to its frozen value ^Tj and it does 
not change a lot between D and B mesons. The reason is that in both cases one has a light 
quark moving in the field of a very heavy one, so that the one-gluon exchange depends on 
the size of the system, rather than on its total mass. Once the difference in size between D 
and B mesons is not that large, the difference between the two values of the strong coupling 
constant is also small (see Table [TVI) . 

The ip function at the origin given in the last column of Table |rV| and which will be used 
later on for spin-spin splittings is defined for radially excited states as 
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where 



(l + A(x-2)), (68) 



POO 

(r^) = (2/ia)^/3^x^) = (2/ia)^/3 / |^^(a;)|2 ^^^^ iV > -3 - 2/, (69) 





which immediately follows from the properties of equation (^5|) and the corresponding re- 
definition of variables. 

B. Spin-spin and spin-orbit splittings. The string correction 

In this subsection we calculate the spin-dependent corrections to the results given in 
Table |IV| as well as those due to the proper string dynamics and which were intensively 
discussed before. 

The eigenstates of the Hamiltonian Hq from (0), which we consider to be the zeroth 
approximation, can be specified in the form of terms ri^^^^Lj (n being the radial quantum 
number) as the angular momentum L, the total spin 5, and the total momentum J = L + S 
are separately conserved by Ho. The corresponding matrix elements for various operators 
present in ( ^5|) read as follows: 



CP^\S,L\^P,) = 0, (iPi|S^2L|iPi) = 0, {^P,\iS,n){S2n)\'Pi) = -|, 

{'Po\S,L\'Po) = -l, ePo\S2L\'Po) = -1, {'Po\{S,n){S2fi)\'Po) = -l, 

(3Pi|5iL|3Pi) = -1, (3Pi|S^2^|3p,) = -1, (^3p^\^s,n){S2fi)\^P^) = l, 

{'P2\S,L\'P2) = \ {'P2\S2L\'P2) = i ('P2\{Sin){S2n)fP2) = ^; 



(70) 
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{'D2\SiL\^D2) = 0, {'D2\S2L\'D2) = 0, {'D2\{Sin){S2n)\^D2) 



1 

"3' 



{^D^\S,L\^D^) = {'D,\S2LfD,) = -^, {^D^\{S,n){S2n)\■'D,) = 



3 /3, 



'D2\SiL\'D2) = {'D2\S2L\'D2) = -I, {'D2\{Srn){S2n)\^D2) = I, 



(71) 



fDal^iLl^Da) = 1, {^D;\S2L\^D3) = 1, {^D3\iS,n)iS2n)\''D,) = ^. 

The interaction Vsd given by (|55|) mixes orbitally excited states with different spins, so 
that the transition matrix elements are given by 



[W2\S^L\^D2) = M {'D2\S2L\^D2) = -M 



(72) 



which lead to mixing within |^Pi), I'^-Pi) and \^D2), |'^-D2) pairs so that the physical states 
are subject to matrix equations of the following type: 



0. 



(73) 



El — E Vi2 
^1*2 E2 — E 

Another important ingredient is the string correction given by ( |52D which leads to an 
extra negative shift for orbitally excited states 



SMi 



-/(/ + l)(r-^). 



(74) 



Gfxjfxl 

Thus the model is totally fixed and the only remaining fitting parameter is the overall 
spectrum shift Cq which finally takes the following values: 

Co{D) = 212MeV, Co{Ds) = UAMeV, Co{B) = 203MeV, Co{Bs) = UAMeV. (75) 

Note that Cq does not depend on the heavy quark (Co(-D) ~ Co{B), Co(-Ds) ~ Co(-Bs)) 
and is completely defined by the properties of the light one. For states with two light quarks 



25 



one would have an overall negative shift 2Co, which gives the contribution — 4Co(M„ — Co) 
to Ml. In case of the p meson this provides a negative constant of order IGeV, i.e., the 
right value needed to bring the theoretical intercept Iq into the correct experimental one (see 
|AM2| in Table imi). 

C. Comparison with the experimental and lattice data. "Mystery" of the D(2637) 

state 

In Tables and we compare the results of our numerical calculations for the spectrum 
and splittings with experimental and recent lattice data as well as with the theoretical 
predictions from and |Q . The underlined figures in Table are considered as the most 



probable candidates for the experimentally observed values. One can see good agreement 
between our theoretical predictions and the experimental values, as well as with the lattice 
calculations p8| , p9| . To demonstrate the relevance of the corrections due to the heavy mass 
we consider a simplified system containing one infinitely heavy particle and the light one 
having its real mass. The best fits for the experimental spectra with the results for such 
simplified systems are also given in Table in the column entitled Mj^^ for comparison. 
One can easily see that corrections in the inverse powers of the heavy mass are strongly 
needed to reproduce the experimental spectrum with a reasonable accuracy and to remove 
degeneracy of S states. 

Now we are in the position to resolve the "mystery" of the D(2637) state (and a similar 
one in the i?-meson spectrum). This state was claimed recently by the DELPHI Collab- 
oration but once its quantum numbers were not defined, then there was a problem 
of the identification of this state. In most quark models (see also Table 0) the first radial 
excitation = 0~ lies approximately in the desired region of mass, but estimates of the 
width of such a state lead to a confusion, as all such estimates give values much larger than 
the width of about 15MeV reported by DELPHI. The only would-be way out of the problem 
is to identify this narrow state with orbital excitations with being 2" or 3^. In spite 



26 



of the fact that orbitally excited states are really narrower than the radially excited ones 
and can have a width compatible with the experimental value, the following two objections 



can be made |41[]: i) quark models predict orbitally excited mesons to be at least bOMeV 
heavier than needed, and ii) a neighboring slightly more massive state should be observed 
as well. 

It follows from Table that we can remove both objections mentioned above (see also 
1^). Indeed, one can easily see that orbitally excited states 2~ and 3~ lie even somewhat 
lower than the radial excitation 0~. The reason for that is the negative string correction 
(ff^) for the orbitally excited states which comes from the proper dynamics of the string. 
Besides that, the single D-wave 3~ state is an even more probable candidate for the role of 
the observed D(2637) resonance than the lightest one from the pair of 2~ states, so that the 
problem of the "missing state" is also avoided. 

Note that our predictions D(2654), Z)(2663), and D(2664) give larger masses compared to 
the experimental one. This must be a reflection of the general lack of the variational einbein 
field method (/i technique) discussed before, which gives slightly overestimated values for 



the spectrum of excited states (see Table | and the discussion in subsection [111 Dj) . 



D. A bridge to the Heavy Quark Effective Theory 

In this subsection we discuss the correspondence between our model and the Heavy 
Quark Effective Theory (HQET) approach widely discussed in the literature (see |]42|j43[| 
and references therein). We use the standard parametrization for the heavy-light meson 
mass, 



Mhi=mQ + A--^{\, + dH\2)+ol^], (76) 



where mq is the mass of the heavy constituent, and the coefficient du describes the hyperfine 
splitting. 
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d 



H 



+3, for states, 
— 1, for 1~ states, 



(77) 



whereas A, Ai, and A2 are free parameters which are subject to theoretical investigation. 
The parameter A2 is directly connected to the splitting between 1^5*1 and 1^5*0 states and 
can be estimated from the experimental 5-mesonic spectrum to be 



A2 = -{Mb* - Mb) ^ 0.12GeV\ 



(78) 



From Table one can easily find our prediction for A2 

A2 ^ 0.16GeV\ 



(79) 



which being slightly overestimated is still in reasonable agreement with the experimental 
value (^. 

In the meantime our model allows direct calculation of the parameters A, Ai, and A2 
based on the Hamiltonian (p^. We apply the variational procedure described above to the 
idealized system with mi = mq —>■ cxo and m2 0. This yields 



A 



-a(Ao) + 



An 



\ 12k 



a(Ao) + 2Ao 



da 






)1 



where k, = |as, and a(A) is the dimensionless eigenvalue introduced in 



Vll A| ); Ao is a solution to the equation 



A2 = f ^ 2A 



da 



dX 



(80) 

(see subsection 
(81) 



Then one can extract the coefficient A2 from the first term in equation (|55D 

Ann 



X, 



3^2 



|^(0)|2 = __o-«:(l + Ao(x-2)). 



(82) 



An analytical formula for Ai is also available, but it is rather bulky and we do not give it 
here. For as = 0.39 one can find the numerical solution to equation ( |81| ) to be Aq = 1.175. 
The corresponding values for A, Ai, and A2 are given in the first column of Table |V11| where 
they are compared with the results of other approaches. 
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Another way to estimate the discussed constants is to find the best fit of the form 

Mj,t = mQ + A + Co--^, (83) 



with A and Ai being the fitting parameters and Co = 203MeV taken from ([75D, for eigenval- 
ues of the Hamiltonian ( |MD with mq varied around the bottom quark mass mi, = A.SGeV 
(we use the region 4GeV < tuq < 6GeV). The coefficient A2 can be found using formula 



with Ao changed for the exact solution for A taken from Table |3 Results are listed in 



the second column of Table Vll 



One can see our figures to be in general agreement with those found in other approaches 
among which we mention the QCD sum rules method [Q, the inclusive semileptonic B- 
meson decays and the Dyson- Schwinger equation for the system of a light quark and a 
static antiquark |^^. We find the parameter Ai to be rather sensitive to the strong coupling 
constant a^. For example, for as = 0.3 one has Ai = — 0.38GeV^^ which should be confronted 
with the value Ai = —0.506GeV^ from the first column of Table |V11| found for = 0.39. 

All our predictions for A2 exceed the value given by equation (0). The reason is the 
slightly overestimated value of ipiO) given by the variational einbein field approach. 

One should appreciate the advantage of the einbein field method which allows one to 
obtain relatively simple analytical formulae for various parameters and to investigate their 
dependence on the strong coupling constant (the dependence on the only dimensional 
parameter a is uniquely restored, giving A ~ y/a and Ai ~ A2 ~ o"). 



VIII. CONCLUSIONS 

In conclusion let us briefiy recall the main results obtained in the present paper. 

We use the model for the QCD string with quarks at the ends to calculate the spectra of 
light-light and heavy-light mesons. There are two main points in which we differ from other 
approaches to the same problem based on various relativistic Hamiltonians and equations 
with local potentials. The first point is that we do not introduce the constituent mass by 
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hand. On the contrary, starting from the current mass we naturally arrive at the effective 
quark masses which appear due to the interaction. Moreover, the resulting effective mass 
is large enough even for the lightest quarks and lowest states in the spectrum, so that the 
spin-dependent terms in the interquark interaction can be treated as perturbations in most 
cases (except for pions and kaons) and thus accounted for in this way. 

The second advantage of the method is that the dynamics of the QCD string naturally 
enters the game and it can be studied systematically. A proper account of this dynamics 
allows one to resolve several problems of the mesonic spectroscopy; namely, one can find 
that the rotating string lowers the masses of orbitally excited states, bringing the (inverse) 
slope of Regge trajectories to their correct values {27ra for light-hght and na for heavy- 
light states). In addition, this allows one to resolve the problem of the identification of 
the D(2637) state recently claimed by the DELPHI Collaboration and which is known to 
lead to a contradiction between its small width, incompatible with the decay modes for the 
radial excitation, and its mass lying considerably lower than the values predicted by the 
quark models for orbitally excited states. In the meantime, taking into account the negative 
string correction contributing into the masses of orbitally excited mesons readily resolves 
this "mystery" for the D(2637) state as well as for its counterpart in the spectrum of B 
mesons. 

For the heavy-light system we extract the constants A, Ai, and A2 used in the framework 
of the Heavy Quark Effective Theory, for which we derive analytical formulae. We find 
our numerical results to be in agreement with those obtained from the experimental data 
and calculations in other approaches like QCD sum rules, inclusive semileptonic S-meson 
decays, and the relativistic Dyson- Schwinger equation for the qQ system. 

We also conclude that the string-like interaction favoured by QCD invalidates the very 
notion of any local interquark potential, still leaving room to the einbein field method for 
the Hamiltonian approach to the bound states of quarks and gluons in QCD. Being rather 
accurate, this method still requires improvements to increase the accuracy and to have the 
full control over it. 

30 



The authors are grateful to A.M.Badahan, P.Bicudo, R.N.Faustov, A.B.Kaidalov, 
V.S.Popov and E.Ribeiro for useful discussions and valuable comments and to V.L.Morgunov 
and B.L.G.Bakker for providing numerical solutions to some equations. One of the authors 
(A.N.) would hke to thank the staff of the Centro de Fisica das Interacgoes Pundamentais 
(CPIP-IST) for cordial hospitality during his stay in Lisbon. 

Financial support of RFFl grants 00-02-17836 and 00-15-96786, INTAS-RFFI grant IR- 
97-232 and INTAS CALL 2000, project # 110 is gratefully acknowledged. One of the authors 
(A.N.) is also supported via RFFI grant 01-02-06273. 



[1] J.Merlin and J.Paton, J.Phys. Gil (1985) 439 

[2] A.Yu.Dubin, A.B.Kaidalov, Yu.A.Simonov, Phys.Lett. B323 (1994) 41; Phys.Lett. B343 
(1995) 310 

[3] Yu.S.Kalashnikova and A.V.Nefediev, Phys.Lett. B492 (2000) 91 
[4] DELPHI Collaboration, P.Abreu et.al, Phys.Lett. B426 (1998) 231 

[5] M.Feindt, O.Podorin (DELPHI Collaboration), contr.paOl-021 ICHEP'96 Warshaw, DELPHI 
96-93 CONF 22; 

C.Weiser in Proceedings of the 28th International Conference on Hight Energy Physics, edited 
by Z.Ajduk, A.K.Worblewski, (World Scientific, Singapore 1997), p. 531 

[6] Yu.S.Kalashnikova and A.V.Nefediev, Yad.Fiz. 60 (1997) 1529 (Phys.Atom.Nucl. 60 (1997) 
1389) 

[7] L.Brink, P.Di Vecchia, P.Howe, Nucl.Phys. B118 (1977) 76 

[8] E.S.Fradkin and D.M.Gitman, Phys.Rev. D44 (1991) 3230; 

G.V.Grygoryan and R.P.Grigoryan, Phys.Atom.Nucl. 53 (1991) 1737; 
Yu.S.Kalashnikova and A.V.Nefediev, Phys.Atom.Nucl. 62 (1999) 377 

31 



[9] A.M.Polyakov, Gauge Fields and Strings (Harwood Academic, Chur, Switzerland, 1987) 

[10] F.Rohrlich, Ann.Phys. 117 (1979) 292; 

M.J.King, F.Rohrlich, Ann.Phys. 130 (1980) 350 

[11] P.A.M.Dirac, Lectures on Quantum Mechanics (Belter Graduate School of Science, Yeshiva 
University, New York, 1964) 

[12] Yu.S.Kalashnikova and A.V.Nefediev, Phys.Lett. B399 (1997) 274; Yad.Fiz. 61 (1998) 871 
(Phys.Atom.Nucl. 61 (1998) 785) 

[13] Yu.A.Simonov, Phys.Lett. B226 (1989) 151 

[14] Yu.A.Simonov, in Proceedings of Hadron'93, Como, edited by T.Bressani, A.Felicielo, G.Pre- 
parata and P.G.Ratcliffe (Nuovo Cim. A107 (1994) 2629); 

Yu.S.Kalashnikova and Yu.B.Yufryakov, Phys. Lett. B359 (1995) 175; Yad.Fiz. 60 (1997) 374 
(Phys.At.Nucl. 60 (1997) 307) 

[15] Yu.A.Simonov and A.B.Kaidalov, Phys.Lett. B477 (2000) 163; Phys.At.Nucl. 63 (2000) 1428 

[16] V.D.Mur, B.M.Karnakov and V.S.Popov (in preparation) 

[17] V.D.Mur, V.S.Popov, Yu.A.Simonov and V.P.Yurov, Zh.Eksp.Teor.Fiz. 105 (1994) 3 (JETP 
78 (1994) 1) 

[18] M.S.Marinov and V.S.Popov, Zh.Eksp.Teor.Fiz. 67 (1974) 1250; 

V.L.Eletsky, V.D.Mur, P.S.Popov and D.N.Voskresensky, Zh.Eksp.Teor.Fiz. 76 (1979) 431 
(JETP 49 (1979) 232); Phys.Lett. B80 (1978) 68 

[19] Yu.A.Simonov, Yad. Fiz. 60 (1997) 2252 (Phys.At.Nucl. 60 (1997) 2069); 
Yu.A.Simonov, J.A.Tjon, Phys.Rev. D62 (2000) 094511 

[20] V.L.Morgunov, A.V.Nefediev and Yu.A.Simonov, Phys.Lett. B459 (1999) 653 

[21] V.L. Morgunov, V.I. Shevchenko, Yu.A. Simonov, Yad.Fiz. 61 (1998) 739 (Phys.Atom.Nucl. 



32 



61 (1998) 664) 

[22] G.F.Chew, Rev.Mod.Phys. 34 (1962) 394; 

I.Yu.Kobzarev, B.V. Martemyanov, M.G.Schepkin, Usp.Fiz.Nauk 162 (1992) 1 (Sov.Phys.Usp. 
162 (1992) 1); 

M.G.Olsson Phys.Rev.D55 (1997) 5479; 

B.M.Barbashov and V.V.Nesterenko, Relativistic String Model in Hadron Physics (Energoat- 
omizdat, Moscow, 1987) 

[23] T.J.Allen, M.G.Olsson, S.Veseli, Phys.Rev. D62 (2000) 094021 

[24] Yu.A.Simonov, Nucl.Phys. B307 (1988) 512; 

Yu.A.Simonov and J.A.Tjon, Ann.Phys. (NY) 228 (1993) 1 

[25] H.G.Dosch, Phys.Lett. B190 (1987) 177; 

H.G.Dosch and Yu.A.Simonov, Phys.Lett. B205 (1988) 339; 
Yu.A.Simonov, Nucl.Phys. B307 (1989) 67 

[26] Yu.A.Simonov, in Proceedings of the XVII International School of Physics "QCD: Perturbative 
or Nonperturhative," Lisbon, 1999, edited by L.S.Ferreira, P.Nogueira and J.LSilva-Marcos 
(World Scientific 2000), p. 60 

[27] A.M.Badalian, Yu.A.Simonov, Yad.Fiz. 59 (1996) 2247 (Phys.At.Nucl. 59 (1996) 2164) 

[28] E.Eichten and F.L.Feinberg, Phys.Rev. D23 (1981) 2724; 
D.Gromes, Z.Phys. C26 (1984) 401 

[29] A.M.Badalian and B.L.G.Bakker in preparation 

[30] J.Pantaleone, S.-H.Henry Tye, Y.J.Ng, Phys.Rev. D33 (1986) 777 

[31] Yu.A.Simonov, Nucl.Phys. B592 (2000) 350 

[32] S.Godfrey, N.Isgur, Phys.Rev. D32 (1985) 189; 
S.Godfrey, R.Kokoski, Phys.Rev. D43 (1991) 1679 



33 



[33] Dan La Course and M.G.Olsson, Phys.Rev. D39 (1989) 1751; 
C.Olsson and M.G.Olsson, MAP/PH/76/(1993) 

[34] A.Le Yaouanc, L.Oliver, O.Pene and J.-C.Raynal, Phys. Rev. D29, 1233 (1984); 

A.Le Yaouanc, L.Oliver, S.Ono, O.Pene and J.-C.Raynal, Phys.Rev. D31 137 (1985); 
P.Bicudo and J.E.Ribeiro, Phys.Rev. D42, 1611 (1990); 
LBars and M.B.Green, Phys.Rev. D17, 537 (1978) 

[35] Yu.S.Kalashnikova, A.V.Nefediev and A.V.Volodin, Yad.Fiz. 63 (2000) 1710 (Phys.At.Nucl. 
63 (2000) 1623); 

Yu.S.Kalashnikova and A.V.Nefediev, Phys.Lett. B487 (2000) 371 
[36] Yu.A.Simonov, Z.Phys. C53 (1992) 419 

[37] Yu.A.Simonov, Yad.Fiz. 58 (1995) 113 (Phys.At.Nucl. 58 (1995) 107); Pis'ma Zh.Eksp. 
Teor.Fiz. 57 (1993) 513 (JETP Lett. 57 (1993) 525) 

[38] J.Hein et.al, Phys.Rev. D62 (2000) 074503 

[39] R.Lewis, R.M.Woloshyn, Phys.Rev. D62 (2000) 114507 

[40] D.Ebert, V.O.Galkin, R.N.Faustov, Phys.Rev. D57 (1998) 5663 

[41] D.Melikhov and O.Pcnc, Phys.Lett. B446 (1999) 336; 
P.R.Page, Phys.Rev. D60 (1999) 057501 

[42] N.Isgur, M.Wise, Phys.Lett. B234 (1989) 113; B237 (1990) 527 

[43] T.Mannel, Effective Theory for Heavy Quarks, in Lecture Notes in Physics, Vol. 479 (Springer, 
New York, 1996), p. 387 

[44] P.Bali and V.Braun, Phys.Rev. D49 (1994) 2472 

[45] M.Gremm, A.Kapustin, Z.Ligeti, M.B.Wise, Phys.Lett. B377 (1996) 20 
[46] Yu.A.Simonov, J.A.Tjon, Phys.Rev. D62 (2000) 014501 



34 



n 





1 


2 


3 


4 


5 


M„ fWKB) 


1.373 


2.097 


2.629 


3.070 


3.455 


3.802 


M„(einbein) 


1.483 


2.256 


2.826 


3.300 


3.713 


4.085 


M„ (combined) 


1.475 


2.254 


2.825 


3.299 


3.713 


4.085 


M„ (exact) 


1.412 


2.106 


2.634 


3.073 


3.457 


3.803 


Af„(combined) — M„(exact) (y 
M„ (combined) ' ''^ 


4.27 


6.57 


6.76 


6.85 


6.89 


6.90 



TABLE I. Comparison of the numerical results of the three approximate methods of solving the 
eigenvalues problem for equation (11) for m = 0, o" = 0.2 GeV^ and / = given by equations 
(p^)-(p5D with the exact eigenenergies of the Hamiltonian Hi from (22). 



n I 


1 


2 


3 


4 


5 





1.719 


2.029 


2.287 


2.516 


2.725 


1 


2.362 


2.611 


2.829 


3.025 


3.208 


2 


2.850 


3.069 


3.264 


3.442 


3.608 


3 


3.259 


3.460 


3.639 


3.803 


3.955 


4 


3.619 


3.806 


3.973 


4.127 


4.268 


5 


3.944 


4.120 


4.276 


4.423 


4.554 



TABLE II. Quasiclassical spectrum of the Hamiltonian (|6l|), ( |62|) for m = and o" = 0.17 GeV^ 
minimized with respect to the einbein /xq (combined method). 
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Meson 






Meccp, MeV 


Mtheor, MeV 


Mtheor [32|, MeV 


Error, % 


p 




1— 


770 


775 


770 


0.6 






3— 


1690 


1688 


1680 


0.1 


P5 


'G, 


5— 


2330 


2250 


2300 


3.4 


Oi 




1++ 


1260 


1346 


1240 


6.8 


as 




3++ 


2070 


2021 


2050 


2.4 


a2 


'P2 


2++ 


1320 


1316 


1310 


0.5 


04 


'F, 


4++ 


2040 


2002 


2010 


1.9 


06 


'He 


6++ 


2450 


2491 




1.7 



TABLE III. Comparison of the masses of the hght-hght mesons lying on the lowest Regge 
trajectories (n = 0, 5 = 1, J = / + 1 for the p and the 02 trajectories; n = 0, = 1, J = / for the 
ai trajectory), calculated for a = 0.17GeV'^ and the overall negative shifts ^/\AM^ = 1126MeV 
for the p trajectory, ^/\AAP\ = 1070Mey for the ai trajectory, and ^/\AM^\ = llOSMeT^ for the 
a2 trajectory, with the experimental data and with the theoretical predictions taken from [^2|. See 



also the discussion concerning the pion and the /o-meson masses in subsection |VI C 
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11 
I V 


7 

I 


TY1 fiQOTl 


TYl -1 
/ / fc]^ 


1112 


U 


U.S 


\ 

/A 


111 

pi 


1 1 r\ 

^2 


1 / 

H- 












B 


1.4 


0.009 


0.17 


0.4 


0.817 


1.497 


0.529 


0.391 


2.198 


0.161 






Ds 


1.4 


0.17 


0.17 


0.4 


0.847 


1.501 


0.569 


0.412 


2.224 


0.167 






B 


4.8 


0.005 


0.17 


0.39 


0.999 


4.840 


0.619 


0.549 


5.527 


0.209 






Bs 


4.8 


0.17 


0.17 


0.39 


1.035 


4.842 


0.658 


0.579 


5.550 


0.219 





1 


B 


1.4 


0.009 


0.17 


0.4 


0.869 


1.522 


0.597 


0.428 


2.640 









Ds 


1.4 


0.17 


0.17 


0.4 


0.891 


1.525 


0.629 


0.445 


2.663 









B 


4.8 


0.005 


0.17 


0.39 


1.052 


4.847 


0.675 


0.593 


5.949 









Bs 


4.8 


0.17 


0.17 


0.39 


1.080 


4.849 


0.707 


0.617 


5.970 








2 


B 


1.4 


0.009 


0.17 


0.4 


0.924 


1.554 


0.674 


0.470 


2.961 









Ds 


1.4 


0.17 


0.17 


0.4 


0.942 


1.557 


0.702 


0.484 


2.982 









B 


4.8 


0.005 


0.17 


0.39 


1.128 


4.860 


0.762 


0.659 


6.245 









Bs 


4.8 


0.17 


0.17 


0.39 


1.151 


4.861 


0.789 


0.679 


6.263 





1 





D 


1.4 


0.009 


0.17 


0.4 


0.929 


1.557 


0.682 


0.474 


2.848 


0.162 






Ds 


1.4 


0.17 


0.17 


0.4 


0.947 


1.561 


0.710 


0.488 


2.869 


0.165 






B 


4.8 


0.005 


0.17 


0.39 


1.142 


4.863 


0.779 


0.671 


6.131 


0.207 






Bs 


4.8 


0.17 


0.17 


0.39 


1.165 


4.864 


0.806 


0.692 


6.149 


0.212 



TABLE IV. Solutions of equations (|6^)-(|67|) for standard values of the string tension a, the 
strong coupling constant a^, and the current masses of the quarks. Eq is the mass of the corre- 
sponding state. All parameters are given in GeV to the appropriate powers. 
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JP 


^exp 


Mtheor 




Mtheor 132] 


Mtheor |4C] 


Mlat [38| 


Mlat H 




± ijy 




1869 


1876 


2000 


1880 


1875 


1884 


1857 


D* 


i _|_ 


\- 




2022 


9000 


2040 


2009 


1 994 


1974 


Di 




1 + 
1^ 


2420 


2354 
2403 


2393 
2407 


2440/2490 


2414/2501 




2405/2414 




l^Po 


0+ 




2280 


2400 


2400 


2438 




2444 


D2 




2+ 


2460 


9439 


9400 


9500 


945Q 




2445 


D*' 




3- 


ZDo ( 


2654 


2600 


2830 










2" 


2663 
2729 


2573 
2600 






0- 


2664 


2500 


2640 


2629 


n 






1968 


1990 


2100 


1980 


1981 


1984 


iTTnn t 

l-Ll Li u 


D* 




1- 


2112 


2137 


2100 


2130 


2111 


2087 


2072 


' - LS 


lip /3 p 


1+ 


ZOoD 


2471 
2516 


2494 
2506 


zooU/zo/U 


zoio/zoOy 


O/IO/I 
Z4y4 


ZOUU/ Zoil 






0+ 




2395 


2500 


2480 


2508 




2499 




13^2 


2+ 


2573 


2547 


2500 


2590 


2560 


2411 


2554 


B 




0^ 


5279 


5277 


5200 


5310 


5285 


5293 


5277 


B* 


± AJl 


1- 


5325 


5340 


5200 


5370 


5324 


5322 


5302 


-'-'1 






5732 


5685 
5719 


5592 
5608 




5719/5757 




5684/5730 






0+ 




5655 


5600 


5760 


5738 




5754 


Bo 


2 


2+ 


5731 


5820 


5600 


5800 


5733 




5770 


B*' 




3- 


5860 


5955 


5800 


6110 










2- 


5953 
6018 


5773 
5827 





38 





2^ Si 


0" 




5940 


5700 


5930 


5898 


5890 




Bs 






5369 


5377 


5400 


5390 


5375 


5383 


input 


Bt 


I'Si 


1- 


5416 


5442 


5400 


5450 


5412 


5401 


5395 


Bis 




1+ 


5853 


5789 
5819 


5793 
5807 


5860/5860 


5831/5859 


5783 


5794/5818 






0+ 




5757 


5800 


5830 


5841 




5820 


B2s 




2+ 




5834 


5800 


5880 


5844 


5848 


5847 



TABLE V. Masses of the D, Dg, B, and Bs mesons in MeV. For the lattice results we give only 



the central values extracted from Figures 26, 27 and Tables XXVIII, XXIX of [38| and from Table 
VIII of ||3^. We also compare out results with theoretical predictions taken from [32| and |4C|. 
The symbols l^Pi/'^Pi and 1^D2/"^D2 are used to indicate that the physical states are mixtures 
of the l^Pi and l^Pi or 1^D2 and l'^D2 states, respectively. Underlined figures give masses of the 
most probable candidates for the experimentally observed resonances. The column Mj^^ contains 
the best fit to the experimental spectrum for the system containing one particle being infinitely 
heavy. 



Splitting 


Ds-D 


D*-D* 


D*-D 


D*s-Ds 


Bs-B 


B*-B* 


B*-B 


B*s-Bs 


Experiment 


99 


102 


141 


144 


90 


91 


46 


47 


Theory 


114 


115 


146 


147 


100 


102 


63 


65 


Theory |32l 


100 


90 


160 


150 


80 


80 


60 


60 


Lattice @ 


100 


92 


110 


103 


90 


90 


30 


29 


Lattice || 


112 


98 


117 


103 


92 


93 


25 


29 



TABLE VI. Splittings for the D, Ds, B, and Bs mesons in MeV. Lattice results are taken from 
Tables XXVIII, XXIX of [|| and Table VIII of pi. We give also results of the theoretical papers 
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mi — > 00 771,2 — ^ 


B mesons 


Sum rules [44 


B mesons decays [45 


DS equation [P[p| 


A, GeV 


0.471 


0.485 


0.4 ^ 0.5 


0.39 ± 0.11 


0.493/0.288 


Ai, GeV^ 


-0.506 


-0.379 


-0.52 ± 0.12 


-0.19 ± 0.10 




A2, GeV^ 


0.21 


0.17 


0.12 


0.12 





TABLE VII. Standard parameters used in HQET (see equation (|76|)). In the first column we 
give the values fohowing from formulae ( ]80| ) and (^) for A and A2 and the corresponding one for 
Ai; the second column contains the best fit of the form iriQ + A + Co — for the eigenvalues 
of the Hamiltonian ( [5^ ) for mg in the region AGeV < mq < 6GeV (around 777;, = 4.8GeV), 
Co = 203MeV being the overall negative shift constant taken from equation (^). The two num- 
bers given in the last column correspond to local and nonlocal kernels of the Dyson-Schwinger 
(DS) equation (see [^^ for details). The figures given in the first two columns correspond to 
a = O.lJGeV^ and a. = 0.39. 



^Note that slightly different values for the string tension (a = O.lSGeV^) and the strong coupling 
constant {as = 0.3) were used in this paper. 
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5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
I.O 
0.5 



U(r) 



1=1 




5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
I.O 



V(r) 



1=2 




lO 



20 25 



30 r 



lO 



20 



25 



30 r 



FIG. 1. Effective potential incorporating the string rotation as well as the quark radial motion 
for a = 0.17 GeV'^ and two values of the angular momentum / (solid line). The naive sum of the 
quark centrufugal barrier for the given / and the linearly rising potential ar is given in each graph 
by the dotted line. For / = the effective potential coincides with ar for all values of r. 




FIG. 2. The Regge trajectories for the Hamiltonian (|56| ) for m = and a = 0.17 GeV^ (see 
Table ||). Theoretical prediction for the /)-meson mass (M^ ~ 1.7 GeV"^; see the column in Table 
I for n = 0) is shown not to violate the straight-line behaviour of the leading theoretical trajectory 
(see left plot). 
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6 r 



5 - 



3 - 



2 - 



1 - 



P(770) 



p^(1690) 



P,(2350) I 



4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



_ «,(1260)I 



a^(2070)J 



5 - 



4 - 



3 - 



2 |_ a,(1320)/ 



%(2040) 



a^(2450) 




FIG. 3. The lowest Regge trajectories for light-light mesons fitted with respect to the overall 
negative mass shift (see Table III). The theoretical values for m = and a = 0.17 GeV"^ are 



marked with dots; the experimental data are given by boxes with error bars. 
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